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Abstract 

We study how the finite-sized n-component model A with periodic boundary 
conditions relaxes near its bulk critical point from an initial nonequilibrium 
state with short-range correlations. Particular attention is paid to the univer- 
sal long-time traces that the initial condition leaves. An approach based on 
renormalization-group improved perturbation theory in 4— e space dimensions 
and a nonperturbative treatment of the q = mode of the fluctuating order- 
parameter field is developed. This leads to a renormalized effective stochastic 
equation for this mode in the background of the other, q ^ modes; we 
explicitly derive it to one-loop order, show that it takes the expected finite- 
size scaling form at the fixed point, and solve it numerically. Our results 
confirm for general n that the amplitude of the magnetization density m{t) 
in the linear relaxation-time regime depends on the initial magnetization in 
the universal fashion originally found in our large-n analysis [J. Stat. Phys. 
73 (1993) 1]. The anomalous short-time power-law increase of m{t) also is 
recovered. For n = 1, our results are in fair agreement with recent Monte 
Carlo simulations by Li, Ritschel, and Zheng [J. Phys. A 27 (1994) L837] for 
the three-dimensional Ising model. 

PACS numbers: 64.60.Ht, 05.70.Jk, 05.70.Ln, 75.40.Gb 
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I. INTRODUCTION 



In the classical work on critical dynamics the important role of initial conditions 
for relaxation processes was not fully appreciated. Guided by the (incorrect) idea that 
all memory of the initial state would be lost after a microscopic time span, so that the 
initial condition would leave no universal traces on macroscopic time scales, most authors 
refrained from a thorough analysis of the effects of the initial condition, taking the relaxation 
process to start in the infinite past t = — oo. Others chose initial conditions at time t = 
and incorporated the initial value rrii = m(0) of the order-parameter density m{t) in the 
scaling analysis of the relaxation process, but tacitly presumed that scaled with the same 
dimension as m(t > 0) @,||J3l- 

As shown some time ago by Janssen et al. this assumption is in general incorrect. 
Due to its failure, quantities relaxing from an initial nonequilibrium state with short-range 
correlations may exhibit interesting short-time anomalies on macroscopic time scales. To 
be specific, let us consider a d-dimensional Ising system with purely relaxational dynamics, 
which relaxes at the Curie temperature T = from an initial nonequilibrium state prepared 
by a rapid quench from a thermal equilibrium state with T ^ Tc and small m^. We assume 
that d is below the upper critical dimension d* (= 4), but above the lower one, d^ (= 2). The 
following is now well-established [Pf-P!^: (i) The initial value rrii gives rise to a macroscopic 
time scale ti ~ mj^^^\ where z and Xj denote the usual dynamic bulk exponent and the 
scaling dimension of mj, respectively, (ii) The latter is greater than = jS/u, the scaling 
dimension of the equilibrium order parameter m(oo), so that 

^,^x^-x^ (1.1) 

is positive, (iii) For t ^ ti, m{t) displays an anomalous increase of the form ~ t^' , which 
crosses over to the familiar nonlinear relaxational behavior ~ t~^^'^^ for t^tj. 



Similar short-time anomalies have been found for other dynamic universality classes |I4 
and in the case of relaxation near a tricritical point [|T^ . 

On a qualitative level, the puzzling fact that m{t) actually increases in the above short- 
time regime even though the equilibrium state that is approached a.s t —>■ oo has zero 
magnetization can be understood as follows (cf. 0). By choice of the initial state, its 
correlation length ^{t = 0) is small. As t grows, ^(t) will increase, too. Fluctuation effects will 
be suppressed as long as ^ (t) is small. Let us suppose for a moment that fluctuation effects 
can be neglected completely for small t. Then a mean-field (MF) description should be valid. 
For systems with short-range interactions this is known to give a critical temperature 
significantly higher than the exact Tc. Hence for short times the system should behave as if it 
were in the ordered MF phase, with m{t) trying to approach the corresponding nonzero value 
of the spontaneous MF magnetization. In reality, fluctuations must not be neglected when 
d < d*, not even in the short-time regime. This is clearly borne out by the renormalization 



group (RG) analyses [@,14], which show that loop corrections do contribute to the short-time 
behavior, producing nontrivial short-time exponents 6' and Xj. Consequently, one must first 
integrate out all degrees of freedom with wave lengths up to the scale ^(t) using the RG 
before reliable information can be extracted by perturbative means. This necessity of taking 
into account fluctuations will lead to quantitative modifications. Yet it seems reasonable 
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to expect that the tendency of m(t) to increase at short times should quahtatively persist. 
Accepting this argument one concludes that as t and ^{t) grow, the system will eventually 
realize that the true is much lower and not below the temperature T of the surrounding 
heat bath. Then m{t) will start its decrease towards zero. 

Recently we have extended the work of Janssen et al. |^ to systems of finite size ||10|,|TT|. 
Specifically, a cube of large, but finite, side length L was considered under periodic boundary 
conditions. As is well known, L gives rise to a characteristic time of the form ti ^ |T^-|Ty[. 
At bulk criticality, this is the only other macroscopic time scale besides tj. In the case of the 
dynamic universality class considered, namely the one represented by the familiar model A 
of Refs. [00], the asymptotic long-time behavior is of the exponentially decaying form 

m(t > ti) ^ Moo e-*/*S (1-2) 

which identifies as the linear relaxation time at T = T^. While is independent of the 
initial state, the amplitude Moo was found to exhibit a universal dependence on both rrii 
and L. It may be written as [|ll| 



Moo(L, m,) ^ L-^/^ g{U/tL) , (1.3) 

where ^ is a universal function, up to redefinitions of its overall amplitude and the scale of 
its argument w = ti/ti. Since Q approaches a constant in the limit w ^ 0, any dependence 
on ti (and hence on rrii) asymptotically drops out completely of (|1.2|) ioi ti ^ ti] i.e., any 
memory of the initial condition is lost for t ^ ^ t^. In the opposite limit, w ^ oo, the 
scaling function behaves as Q{w) ~ ^ooW^"^'^^, giving ||I0[ 



Moo{L, rrii) ^ const x rrii L""^' ■ (1-4) 

That is, the amplitude varies linearly in mj, and its L-dependence is governed by the same 
exponent 6' that characterizes the short-time singularity of m{t) in the bulk case. 

As we have shown in Ref. these predictions follow in a natural fashion from a 



phenomenological scaling approach. While they should therefore be valid for general values 
of n, detailed checks and explicit calculations of the scaling function Q have so far been 
undertaken with only two methods. One, utilized in Ref. |T0[, is the exact solution of the 



finite-sized n-component model A in the limit n —>■ oo. This yielded for the amplitude in 
( |1.2D the simple form 

Mr-^ = , (1.5) 



where A and B are nonuniversal constants. 

The second method used so far is Monte Carlo simulation of two and three-dimensional 



Ising models [T^|T^. Both Li et al. and Grassberger |]T3[ estimated the short-time 
exponent 6', obtaining values in reasonable agreement with the e-expansion results of Janssen 
et al. M. In Ref. |T^, also the dependence of M^o on rrii was investigated and found to 



conform remarkably well with the large-n result ( p..5| ). 



In the present paper we wish to apply an alternative powerful approach, one which 
works for general n: the RG approach that uses e = 4 — d as small parameter. From 
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previous studies P,[7|,p!7|-[T9[] it is well known that a proper treatment of finite-size effects by 
this method requires a rather subtle procedure. Unlike to the case of bulk critical behavior, 
a dimensionality expansion in integer powers of e does not in general result. The reason 
is that the lowest-energy mode of the order-parameter field 0(a;,t), which is the spatially 
homogeneous one (with wave number q = 0) in our case of periodic boundary conditions, 
must be treated nonperturbatively; only the other, q modes may be treated by RG- 
improved perturbation theory. 

A frequently used technique P,p!8|- pl|] is to integrate out the latter modes in order to 
obtain an effective action for the former. Upon using the functional-integral formulation of 
our dynamic model, this could in principle also be done in the present case. Due to the 
fact that even the free propagator of the q ^ modes involves the q = part of as a 
background field, the resulting effective action becomes rather complicated, however, and 
does not lend itself easily to direct calculations of time-dependent quantities such as m(t). 
We have found it more convenient to perform the perturbative elimination of the q ^ 
modes directly on the level of the Langevin equation, following Goldschmidt's procedure 
in 1^ in this regard. At lowest nontrivial order of renormalized perturbation theory the 
resulting equation for the q = mode takes the form of a nonlinear stochastic equation. 
The required nonperturbative treatment of this mode then amounts to solving this equation 
by numerical means. Upon averaging the solutions over the associated noise, the desired 
m(t) follows. 

Our analysis goes in two respects beyond the related one by Goldschmidt 0. First of 
all, we take care to properly build in the initial conditions. This is a central purpose of the 
present work. The advantage is twofold. On the one hand, our approach encompasses an 
adequate description of the anomalous short-time behavior. On the other hand, it permits 
us to study in a systematic fashion how the initial condition rrii manifests itself on long time 
scales. In this regard, our work complements — and in parts corrects — Ref. 0. 

Second, we consider the finite-sized model A for general component number n. Accord- 



ingly our findings can be compared both with the exact results for the limit n ^ oo |10,11 



as well as with Monte Carlo calculations for the kinetic Ising model (n = 1) ||T2[ and for 
polymers {n = 0) p!3| . 

The remainder of this paper is organized as follows. In Sect. II we first recall the definition 
of the finite-sized n-component model A. Starting from the stochastic equation for the 
order-parameter field (f){x,t), we eliminate the q ^ degrees of freedom to obtain the 
reduced stochastic dynamics for the q = component of 0. Section III deals with 
the finite-size scaling forms of the stochastic equation for and of its solution. We 
first reconsider briefly the problem in the case d > 4. Upon renormalizing the theory for 
d < 4, RG equations are then derived for <I>(t) and its stochastic equation whose solutions 
yield the desired asymptotic scaling forms. Subsequently the so-obtained scaling form of 
the stochastic equation is explicitly verified by means of RG-improved perturbation theory 
up to one-loop order within the framework of the resulting y/e expansion. In Sect. IV we 
first show that the scaled stochastic equation at bulk criticality has asymptotic solutions 
with the correct limiting behavior in both the short-time regime as well as in the long-time 
one. The scaled stochastic equation is extrapolated to e = 1, and numerical results for its 



solutions are presented. These are compared with the large-n results of Ref. [10| and the 



Monte Carlo simulations of Refs. and [|T^. Section V contains a brief summary and 
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concluding remarks. Some technical details have been relegated to two appendices. 



II. DEFINITION OF THE MODEL AND REDUCED DYNAMICS FOR THE 

ZERO-MOMENTUM MODE 

We consider the ra-component model A for a d-dimensional hypercube V with volume 
L'^ under periodic boundary conditions. This is defined by the Langevin equation 

9,0(a;,t) = -Ao|?|^ + C(^,t) (2.1) 



with the Ginzburg-Landau Hamiltonian 



n{4>} 



V 



(2.2) 



where = (0") is an n-component field, while ( = is a Gaussian random force with 
mean zero and correlation 

{C{x, t) C^{x', t')) = 2Ao 5"^ 5ix - x') 5it - t') . (2.3) 

All times t and t' are assumed to be positive. At t = we impose the initial condition 

0"(a;,O) = 5°^mi , (2.4) 

where rrii is independent of x. 

In view of the chosen periodic boundary conditions we make a Fourier expansion 



(a.,t) = L-'^$:0^(t)e^''-, (2.5) 



where the sum runs over all q = y^m with m E Next we decompose 



r{x,t) = (i)^{t) + rix,t) (2.6) 

into its g = part 

roit) = L-'' j (2.7) 

and a remainder = (?/)") given by the sum of all g 7^ terms on the right-hand side of 
(p.5|). An analogous decomposition we make for the noise, writing 



ax,t) = E{t)+i^{x,t) . (2.8) 
As a consequence of (p.3|) the q = part 

E{t)^L-' f C(x,t) (2.9) 



is a Gaussian white noise with correlation 
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S^(t')) = 2Ao L-'' 5{t - i!) . (2.10) 
The 9 7^ parts ^/^ and satisfy 

'"Xa;,t) = 0. (2.11) 



Y 



By averaging ( p^.l| ) over the 9 7^ part -(9 of the noise we can formally obtain a stochastic 
equation for the constant mode alone. Let us denote an average over d by (.)^. Further, we 
introduce the notation 

0-^$"(t) + 50"(t), (2.12) 



with {6(f)'^{t))^ = 0. Because of ( p.ll| ) and translational invariance we also have {ilj{x,t)) 



1? 



0. The resulting stochastic equation for $(t) can be written in the form 



<i>"(t) + Aoro$"(t) + Uf^(mvm\t) + 3$^(t) C^\t) + Cf^^^t)] = H"(t) , (2.13) 

6 L -I 

where 



and repeated indices are summed over. A systematic derivation of Eqn. ( |2.13| ) is given in 



Appendix A. The functions C"^(t) and C°'^'^{t) in (|2.14|) are the expectation values of bilinear 



and trilinear products of ijj^s and 6(f)^s. These are independent of x but depend, of course, 
in a complicated manner on $(t') with < t' < t. In terms of $ the initial condition ( p.4|) 
becomes 

$"(0)=r^mi, (2.15) 

and expectation values involving ip^s and 50's have to satisfy Dirichlet initial conditions. 
Eqn. ( p.l4|) , which is an exact consequence of ( |2.1| )- ([27^ ), forms the basis of our sub- 



sequent analysis. We emphasize that the functions C"^ and C"^"^ are functionals of $. In 
order for (|2.1|) to become a closed equation for $, this functional dependence must be known. 
This is the case — at least in principle — if we consider these functions to be given by their 
Feynman graph series. Explicit expressions for their graphs to two-loop order are presented 
in Appendix A. The ^-dependent free propagators involved — and hence the resulting ex- 
pressions for the graphs — differ from the ones used by Goldschmidt [0. This is because 
we imposed an initial condition at t = rather than in the infinite past. 



III. SCALING FORM OF STOCHASTIC EQUATION AND OF ITS SOLUTION 

We are interested in the asymptotic solutions of the stochastic equation ( p.l3| ) near the 
bulk critical point. Since these are known to have a scaling form, the stochastic equation 
itself should take a scaling form on sufficiently large length and time scales. In order to show 
this and to compute explicit expressions for these scaling forms, we will largely follow Ref. 
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Differences arise because of our imposed initial condition at t = 0. In the case d < A, 
essential modifications are required to ensure a proper treatment of both the short-time 
anomalies and the long-time behavior. In the case d > 4, the differences are marginal as far 
as the leading asymptotic behavior is concerned, and the results could be gleaned from Ref. 
0. Since the subsequent analysis for d < 4 makes substantial use of these, it is useful to 
recapitulate them briefly. 



A. The case d > A 

Above the upper critical dimension d* = 4, the effect of the q modes (i.e. of tp) 
is asymptotically negligible. Thus C"^ and C"^'^ may be replaced by zero, their zero-loop 
expression, in ( |2.13| ). Inspection of the resulting stochastic equation reveals that it takes a 
scaling form. In summary, it follows that $ asymptotically behaves as 



m^%'^'L-^^'M ^,^,^o/V/V;^^M , (3.1) 



where the scaling function $ = ($") is a solution to the scaled stochastic equation 

1 



at + l + -|$(t)|^ 



$"(t) = S"(t) (3.2) 



with the initial condition 

<|"(t = 0) = b''^g]I^L'^l^m^ . (3.3) 

Here t = gl^'^L^'^^'^Xot is the scaled time, and the scaled noise H(t) = XQ^gQ^^^L^'^^^E{t) 
obeys 



;S"(t) H^(t')) = S{t - t') . (3.4) 
The result (|3.1| ) is a manifestation of the familiar fact that finite size scaling p^j23[1 gets 



modified above the upper critical dimension [p4| , P5i [4 [r7|[r9| . Below it the finite size scaling 
forms involve the temperature field r ~ (T — Tc)/Tc and the time t in the combinations 
C,/L and t/L^, respectively, where C, ~ t^*^ is the equilibrium bulk correlation length. This 
simple form of finite size scaling breaks down for d > 4 because go gives rise to a further 
length that must be taken into account in the scaling theory. 

A second and trivial remark concerns the scaling of m^. Since the scaling form ( p.l|) is 
dictated by dimensional considerations, necessary scales with L in precisely the same 
fashion as $ and hence m{t). Only for < 4 can the corresponding scaling dimensions Xi 
and Xfj, become different. 

Finally, we note that corrections to the leading asymptotic term given on the right- 
hand side of ( |3.1| ) can be computed by perturbation theory in go, provided the ultraviolet 
singularities have been regularized. 
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B. The case d < 4 



1. Renormalization 



In this case the influence of the q ^ part ip on ^ must not be neglected. Since 
perturbation theory in qq is known to break down in the critical regime, it cannot be used to 
treat this influence. We will rely on RG-improved perturbation theory near d = 4. To this 
end we must first renormalize the theory. Utilizing the scheme of dimensional regularization 
and minimal subtraction of poles in e, we introduce renormalized quantities through 



X, t) = Z'J' t) , (f){x, t) = Z'f t) 



1/2 



(3.5) 



2 -2 / / \ 

Tq- To^c = ^^ Z^T , goKd = iJ^'Z^u , Ao = /i [Z^/Z^j X 



(3.6) 



t=o 



,1/2 



t=0 



R 



and 



C(a;,t) = (V^^)'^'Ci?(a^,t)- 



1/2 



t=0. 



R 



(3.7) 



(3.8) 



Here is the auxiliary field that is needed in the functional-integral formulation of 
the dynamics specified by (|2.1j)-(p7^). We have absorbed the usual angular factor = 
2/[(47r)'^/2r((i/2)] in u. The momentum scale /i is arbitrary and will be set to /z = 1 in the 
sequel whenever it is suppressed. Further, Z^, Z^, Z^-, and Z^ are standard renormalization 
factors. Finally, Zi is the analog of the renormalization factor denoted Zo by Janssen et al. 



For the associated RG functions we use the notation 

ri^{u) = fidf,\^\nZ^ , k = (j), (p, t, u, i 



and 



(3.9) 



(3.10) 



where d^l^ means a derivative at fixed bare variables. Using familiar RG arguments, the 
critical exponents can be expressed in terms of the values rj* of the functions ( |3.10| ) at the 
infrared-stable zero u* of Specifically for the scaling dimensions and Xi and the critical 
exponents z and 9' one has [H 



d-2 + v; 



^1 =i< I ^ 



^1 =f= I ^ 



2z 



(3.11) 



In the one-loop calculation presented below we shall need the 0{u) results 
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l + + Z^ = l + '^u + 0{u') , (3.12) 

be be 



Z^ = l + 0{u'), Z^ = l + 0{u^), (3.13) 



and the e expansions 



u* = -^ + 0{e') (3.14) 
n + 8 



and 



We must also discuss what happens to the initial condition ( |2.15D upon renormalization. 
This condition can be incorporated into the stochastic equation (|2.13| ) by adding the term 
(5(t)mj$^(t) on the right-hand side. In the equivalent functional-integral representation of 
the theory, such a term corresponds to a contribution of the form — /y 0^(a;, 0) to the 
action (cf. Eq. (5.28) of Ref. P]). The fact that the renormalization ( |3.7|) of 

1 /2 

an extra factor implies that nii should be reparametrized as (see Eq. (67) of Ref. 



involves 

t=o 



— 1 /2 

[Z^Z^ nii^R . (3.16) 



rrii 



2. RG equations and finite-size scaling forms 



In this part, we will only be concerned with the renormalized fields and Sr, but not 
with the bare ones. To simplify the notation, we shall drop the subscript R on these as well 
as on mi^R. Thus $, S, and now mean renormalized quantities. 

By exploiting the arbitrariness of /i in the above reparametrization relations, RG equa- 
tions for $ and its stochastic equation can be derived in a standard manner. Let us write 
the (renormalized) stochastic equation for $ as 



/C(t,L,r,u,A,/i,{<l>,S}) = 

and introduce the operator 

'D^, = fidf, + Pudu - (2 + r]r) rdr + ^ (4 + r?^ - r]^)\dx ■ 
Then the RG equations for $ = ($") and /C = (/C") can be written as 



T^a^ ^ "^i7^ 

^2 2 dm, 



$ = 



(3.17) 



(3.18) 



(3.19) 



and 
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2 



/C = 



(3.20) 



We solve these by characteristics. The solutions involve metric factors, given by inte- 
grals of exponentials of functions along RG-trajectories. To account for the nontrivial 



e-dependence of the scales implied by the ^fQ-dependence of the zero-loop result (|3.21| ), it is 
convenient to split off appropriate powers of a constant ~ e from the metric factors. We 
leave open the choice of at this point; the choice ( p.43[ ) we will finally make is such that 
the effective coupling constant appearing in the scaled form of the stochastic equation of 
motion when T = approaches one in the long-time limit. 

The resulting asymptotic finite-size scaling forms then become 



L, m,, r, u, A, {S}) ^ G^'^^ L'^^ <l> (t, t„ y, {h}) 



and 



/C(t, L, r, u, A, /i, {$, H}) ^ XGl/%'L-^/^ }C{t, y, {$, s}) , 
with the familiar magnetic scaling index A/u = d — and the scaled arguments 

t = tXGl'^BxL-' = t/tL , 



(3.21) 



(3.22) 



(3.23) 



t, = (S,m,)-^/^' Gl/'B,L-^^U/t, 



(3.24) 



and 



y = G-'/'BrTL^''^ 



(3.25) 



(3.26) 



where = B^{u), Br, Bx, and Bi are nonuniversal metric factors. Consistency requires 
that the scaled variable $ in ( ^.22| ) be given by 



(3.27) 



For the second moment of the scaled (renormalized) noise S we recover ( |3.4| ). The scaling 
functions $ and /C are given by the functions $ and /C with X = L = fi = 1 and u = u*, i.e., 



/C(t, y, {$, H}) = G;'/'' IC{G:'/\ 1, T = G'J'y, u\ 1, 1, {G;i/^<l, Gl''^]) 



(3.28) 



Below we shall explicitly verify the scaling forms ( |3.21 ) and ( 3.22 ) by means of RG- 
improved perturbation theory to one-loop order. This will be done in two steps. First, we 
shall compute the scaled equation of motion K, for arbitrary y > 0. Then this will be solved 
numerically at the critical temperature to determine the scaling function $ for y = 0. 
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3. Explicit calculation of the scaling form of the stochastic equation 



Our starting point is the stochastic equation ( |2.13D for the bare field $. (Since the bare 
and renormahzed fields $ and S agree at the one-loop order we will be working, there is no 
need to restore the subscript R on the renormalized ones. Whether $ and S mean bare or 
renormahzed fields should be clear from the context anyway.) 

In order to be able to apply RG-improved perturbation theory, we first must express all 
bare quantities in terms of renormalized ones. We begin by substituting the loop expan- 
sions of the functions C"^ and C'^^"' discussed in Appendix A into ( p.l3|) , where we retain 
only contributions up to one-loop order. According to this appendix, in which the explicit 
expressions for the Feynman graphs of C"^ and C"^'*' are given to two-loop order, the latter 
function vanishes at one-loop order and hence can be dropped. The one-loop result for the 
former reads in matrix form 



c{t) = 2X0 r 



dt' 



1 



— 2AorQt' 



(3.29) 



with 



A{x)= J2 e 



(3.30) 



and 



£{t', t) = exp 



■^Ao^o f dt"M{t-t" 

3 JQ 



where the matrix Al(t) is defined by 

Wf^it) = $2(t) 6""^ + 2<l>"(t) $^(t) 



(3.31) 



(3.32) 



Note that the t' integration in ( p.29|) ranges from to t. In the analogous Eq. (3.24) of 
Ref. [0 this integration starts from —00. This is the reason why the anomalous short-time 
behavior was missed. 

A second observation is that the one-loop term ( ^.291 ) contains an ultraviolet singularity. 
To see this one should recall from the analyses 0,^,^] of the static theory that the function 
A{x) varies as A{x) ~ (vr/x)^/^ as x — 0. Thus the integrand in ( p.29|) behaves as t'~2+e/2 
near zero, so that the t' integration diverges at its lower limit. 

In order to extract the associated pole parts and to show that they have the required 
local form, we proceed as follows. Suppose we expand S{t', t) in powers of t'. Then only the 
terms to linear order in t' will contribute to the pole parts. Accordingly we subtract and 
add in ( ^.29] ) a term corresponding to this Taylor expansion to order t'. The analog of C{t) 
of C(t) defined through replacement of S(t',t) by 



e{t', t) = Sit', t) - 1 + -AoTo t'M{t) 



(3.33) 



in ( p.29|) is regular in e. The added terms suggest the definition of the time-dependent (and 
L-dependent) variables 
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ro{t) = To + 



exp (2Aorot') 



and 



n + 8 2Xlgi ft ^, A'^i8n'X,L-H') - 1 



exp (2Aorof ) 



(3.34) 



(3.35) 



In terms of these quantities the stochastic equation that holds for $(t) at one- loop order 
can be written as 



S""^ + ^Ao^o ^"^^^C^^(t) |>$^(t) = H°(t) . (3.36) 



The Laurent expansions of the integrals appearing in ( 3.34| ) and ( p.35 ) are computed in 
Appendix B. The results are given in (|B8D and We substitute these into (|]3|) and 

( p.35| ). Upon expressing go, tq, and Aq by means of ( p.5|) -( p^ and ( p.l2| ) in terms of their 
renormalized analogs, the one- loop poles ~ u/e in fo(t) and go(t)/u are found to cancel, 
giving 



fo(t) = m + Oiu') , goit) = m + Oiu') 
with the ultraviolet-finite, time-dependent quantities 



n + 2 I 



Inr - Ei(-2Art) 



exp(— 2Atr) 



2Art 



(3.37) 



(3.38) 



and 



~g{t)=K^'\u + 
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1 + lnr-Ei(-2At 



r 



+ 0{e) 



(3.39) 



where Ei(x) is the exponential-integral function in the notation of [^, and we have intro- 
duced the functions 



Fp{a,(3) = rdxxP (^A\x) - 1 - n^/x^) e 



-I3x 



(3.40) 



with p = 0, 1. The quantities f and g are time- dependent analogs of the static shifted 
variables introduced in Refs. PO] and encountered again in the treatment of the dynamic 
case in Ref. [0 (where a different normalization of the renormalized coupling constant was 
employed, however). In the limit t —>■ oo, our t-dependent variables f and g reduce to these, 
as they should. That is. 



n + 2 



■ U T 



4 / 

\nr + —Fo[oo,L^T/4n' 



(3.41) 



12 



and 



lim g{t) = g^t = K^^ull 



t^oo 



n + I 
~l2 



u 



1 + lnr r Fi(cx),LV/47r 



(3.42) 



In the limit t — > 0, f(t) and g{t) have short-time singularities of the form f{t) ~ t^^ + 
(9(rlnt) and g{t) ~ Int. As we will see below, the singularity has the correct form and 
coefficient to produce the short-time singularity ~ of 

The renormalized stochastic equation that holds for at one- loop order is given by 
( p.36| ), with To, S'O) and Aq replaced by f, g, u/K^, and A, respectively (and $ and S 
interpreted as renormalized fields). We now choose the constant in (|3.21|) - (|3.28|) as 



= lim gst{r,L = l,u* 

T-+0+ 

n + 8 



KJ\* { 1 + 



12 



u 



l-Cs + ln47r2-i^/i 



(3.43) 



with 



/i 



A{x) 



1 — TT X 



2-2 



[l-x) 



(3.44) 



where Ce is Euler's constant and 6{x) denotes the Heaviside step function. 

Upon substitution of the above results into (|3.28| ) the scaled stochastic equation /C = 
becomes 



5t + reff(2/,t) + ^Ucff(y,t)|<|)p 




$ = S 



(3.45) 



with 



reff(t,y) = G';i/'W 



12 



In (yG'J') -Eii-2yt) 



exp(-2yt) 



2yt 



+ 



n + 2 



3/2 



(3.46) 



and 



1 ~* T z—2x 



Uesit, y) =G^ g* L 



1 + 



n + 8 
12 



u 



GE + \n {Gl/'y /47r')-Eii-2yt) 



O e 



(3.47) 



The contribution involving 0°"^ would give rise to a memory term; this has been omitted 
since it is of order e^^^. 

These equations hold for general values of t and y. However, the numerical calculations 
discussed in the next section have only been carried out at the critical point, ?/ = 0. For this 
value of y the expressions (|3.46|) and ( p.47|) simplify to 
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.5,(t) . f^i^, 0) ^ ^ v^F„(8.^G.-V.,o) - !i±| 1 i + (3.48) 

and 

00 

ulsit) = 7ieff(t, 0) = 1 + ^ / dxx [A{xy - 1] + 0(e2) . (3.49) 

These quantities approach nonvanishing constants f^s{oo) and -1(35(00) = 1 in the hmit 
t ^ 00. In the short-time hmit they behave as 

and 

u^sit ^ 0) ^ 1 - I [in {8n'G;'/h) - tt-^i + 0(t)] + 0{e') , (3.51) 

respectively. As we will see below, the short-time singularity of the latter quantity should 
be exponentiated to a power of the form i-<^/'^+^i^ ). 

IV. NUMERICAL SOLUTION OF SCALED STOCHASTIC EQUATION 

We now turn to the computation of the finite-size scaling form of the order parameter 
m(t). To this end we must solve the scaled stochastic equation (|3.45|) and average the 



solution $(t) over the scaled noise S(t). Restricting ourselves to the case of relaxation 
directly at the critical point, we set r = y = 0. For notational simplicity, we also put the 
time unit A^^ = 1. Further, we fix the nonuniversal metric factors and Bx such that 
G* ^^'^B^ = gV^Bx = 1, leaving open the choice of Bi for the moment. 

Upon averaging the scaling form (|3.21|) of $^(t), we see that at the critical point m{t) 



can asymptotically be written as 

m{t, L, mi) ^ L-"^ a(t, t^) = Bim^ f f{t, U) (4.1) 

with 

a(t,t.) =t^r/7(t,t.) = (l>^(t,t„0,{H}))„ (4.2) 

where (.)a means an average over S. The advantage of rewriting the scaling form in terms of 
/ is that the short- distance singularity has been split off. Hence limt^o /(t, U) must exist 
and be nonvanishing, and since m{t, L,mi)/mi must be independent of whenever t <S 
min(ti, tj), the limiting value /(O, t,) must be independent of t,. We choose the normalization 

/(0,t,) = /(0,0) = l, (4.3) 

which fixes the nonuniversal constant Bi. Upon redefining rrii, we absorb this constant in 
rrii, and hence take ;Bj = 1 in the sequel. 
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In general, the solution of the scaled stochastic equation (p.45|) can only be obtained by 



numerical integration. However, some analytical information about the asymptotic behavior 
in certain regimes can be extracted. Let us consider ( |3.45| ) for small t = t/ti. Due to the 



short-time singularity of f^g, all bounded continuous solutions of ( p.45| ) must start from zero. 
(No confusion should arise from the fact that $^(t, ti) and hence the scaling function o"(t, ti) 
in (|4.1| ) approach zero as t ^ even though mj 7^ 0. In view of the anticipated short-time 
singularity o" ~ and ( |1 . 1[ ) this is seen to be an immediate consequence of the inequality 

Provided rrii was taken small, <l>^(t, tj) will also be small for t <^ min(tj^,tj), i.e., for 
t ^ min(l,ti). Hence the nonlinear term may be neglected in this short-time regime. The 
solution to the resulting linear equation (involving the short-time singularity of is indeed 
of the form $^(t, tj) oc , with 9' given by the 0{e) expression ( p.l5| ). 

As t grows, ^ grows. Since the proportionality constant is proportional to mj, this initial 
growth is very rapid if mj is large. Once t reaches tj, the nonlinear term in the stochastic 
equation ( p.45| ) must no longer be neglected. The larger t/U gets, the more important this 
term becomes till it finally dominates completely both the linear term oc f^g and the noise. 
As long as t remains small in comparison to ti-, we may then use the asymptotic short-time 
form -Ugfj ~ t~^/^+'^('^ \ In this manner one finds that the solution decays as 

a(t< l,t, <t) -t-'^/'^^ (4.4) 

in the specified regime, where the exponent which is known to govern the nonlinear 

critical relaxation of the order parameter in bulk systems, takes the one- loop value (3/vz = 
1/2 - 6/4. 

In order to determine the asymptotic behavior in the long-time regime t ^ 1, one can 



substitute t^^ and in (|3.45|) by their limiting values rgfj(oo) and uIq{oo) = 1. The 



resulting stochastic equation is equivalent to the one investigated in previous analyses of 



finite-size effects on linear relaxation 1^,^,17]. From these papers and Ref. |10] it is clear 
that o"(t, tj) decays as 

a(t^oo,t,)^oroo(t.)e-^°"^*^^ (4.5) 



where the amplitude o"oo(tj) according to Ref. [jlO[ should depend on tj in a nontrivial fashion 
(cf. the discussion in the Introduction). In particular, it should have the asymptotic behavior 

, , f a for tj — > , , , 

[ 6 tj ' = 6 rrii for ti ^ 00 , 

where a and b are nonvanishing universal constants, and 

m. = tr/^ = m,tI^/\ (4.7) 

We expect that it should be possible to extract analytical information about croo(ti) from 
the scaled stochastic equation ( p.45|) . We have not done this, however. In the following we 
concentrate on the numerical solution of (|3.45D, which provides detailed information about 



all regimes discussed above. 

For the numerical solution of (|3.45| ) in = 3 dimensions we make the ansatz 
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(4.8) 



with 

This leads to the equation 
d 



1 



^ + <frM + ^M^ff(t)|$P 



1 . 



with 



6' n + 2 



cffV 



and 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



where denotes the n — 1 dimensional component of $ in the subspace orthogonal to the 
a = 1 direction. 

For reasons discussed above it is important that the short-time singularity of u'^s{t) is 
properly exponentiated. We used the e = 1 expression 



^eff(t) = 1 + 



1 



VOAx 



-0.62: 



X 



(4.13) 



which smoothly interpolates between the asymptotic form u^g ~ t^*^/^ at small t and the 
large-t limit 'Ugg(oo) = 1 (see Fig. 1). The constants 0.4 and 0.6 in ( ^.13| ) were determined so 
that the interpolation formula ( |4.13 ) fits the the numerical results of the expression (|3.51 ) 
with e = 1 in the large-t regime . For r^g(t) we utilized the expression (|4.11| ) with e = 1. 




Fig. 1.: 



0.2 0.4 0.6 0.8 1 



as a function of x 



Svr^G* t. The solid line represents 
the data from the numerical integra- 
tion of ( 3.4S| ) with e = 1. The dashed 
line shows the asymptotic logarith- 
mic behavior for a; — > 0, given in 
( p. 51 ). The dotted curve repre- 
sents the expression (111), which 
was used in the numerical calcula- 
tion. 



The stochastic equation (|4.10| ) was integrated numerically with a first-order Euler pro- 
cedure ||2^. For n ^ 1 one has to solve simultaneously also the original equation (p.45|) 
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for the n — 1 orthogonal modes. Averages for the orthogonal modes remain at zero, which 
means that the order parameter does not 'rotate' within the n-dimensional space. In order 
to obtain averages, summations over 30 000-100 000 single trajectories were performed. Ex- 
emplary results for a{t, tj) with rrij = t^"^"^^ as parameter are displayed in semi-logarithmic 
and double- logarithmic representation in Fig. 2 for n = 1 and in Fig. 3 for n = 20. 




Fig. 2.: Order-parameter profiles for 
n = 1 and rrij = 0.1, 0.2, 0.5, 1, 2, 
5, 12, 25 (from bottom to top) ob- 
tained from the numerical solution of 
(P^ ). In the main figure the data 
are displayed in semi-logarithmic form. 
The small double-logarithmic diagram 
inserted shows the same data empha- 
sizing the power-law growth for small 
time and the nonlinear regime for 
larger values of m, (upper curve). The 
pure power laws ~ t^^^^ and ~ t^^/^ 
are represented by the dashed lines. 




Fig. 3.: Order-parameter profiles for 
n = 20 and ruj = 0.1, 0.2, 0.5, 1, 2, 
5, 12, 25 (from bottom to top). The 
dashed lines in the small figure show 
the pure power laws ~ ^1^/56 gj-^^ ^ 

^-1/4^ 



A common property of all curves is the expected power-law increase for small t. The 
lower dashed lines in Figs. 2 and 3 represent the pure power laws t^/^^ and t^^^^^, which, 
according to ( p.l5| ), hold for n = 1 and n = 20 in the limit t ^ at the level of our 
approximation. In the large-t regime the scaled magnetization decays exponentially in all 
cases, with an rrij-dependent amplitude. For small rrij (bottom curve), the initial algebraic 
growth crosses over directly to the linear relaxational decay. As rrij increases, the maxima 
of the curves gradually shift to the left. Further, immediately to the right of the maxima, a 
regime of nonlinear relaxational behavior starts to emerge, which becomes more and more 
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pronounced as rrij grows. 

As the results in Figs. 2 and 3 show, the scahng function a{t,ti) becomes independent 
of rrij in the hmit rrij — > oo (tj — > 0). Its hmiting form o"(t, 0) has the following properties: 
At t = 0, it jumps from to oo, then it exhibits nonlinear relaxational decay ~ t~^^'^^ until 
the crossover to linear relaxational behavior occurs. This is the special case studied in the 
existing literature [0. The top curves in both figures are already close to this limit. They 
approach the upper dashed lines, which represent the power law t~^^^, 1/4 being the value 
of jS/uz in our approximation. In the opposite limit rrij — (tj —>■ oo), the scaling function 
cr(t, tj) approaches, of course, zero for all t. 




We have also determined the amplitude o"oo(tj) defined in ( [4.5|) . The results are displayed 
in Fig. 4. They are in conformity with the limiting forms for large and small U predicted in 
( [4.6|) . It is interesting to compare them with the exact analytical expression 

^^r^\U = mr^") = , ^7 ^ (4.14) 

[a/b] + 



that holds in the limit n oo [|1^]. To this end we have tried to fit our numerical data 



to the form ( [4. 141) . For n = 20 the function ( [4.14|) provides a good description of the 



numerical results (represented by the squares) . For the curve shown in Fig. 4 the parameter 
values are a = 4.92 and b = 0.96. For n = 1, however, ([4.14|) yields a relatively poor fit, 
especially in the crossover region. The numerical data are represented by the full circles. 
The crossover region is magnified in the small diagram inserted in Fig. 4., where the dashed 
curve shows the fit obtained with ( [4.141 ). A much better fit is provided by the function 



a rrij/ {{a/bY'^ + mf-^)^^^'^, with a = 2.38 and b = 1.13, depicted by the solid curve. 

In our discussion of the short-time anomalies of critical relaxation in the Introduction 
we pointed out that this phenomenon is intimately related to the fact that the correlation 
length ^{t) is initially small and grows as time proceeds, approaching its equilibrium value 
^eq = C,{oo). Rather than the correlation length ^(t), we preferred to investigate the variances 
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m(2)(t) = - m{tf , (4.15) 

and 

m^l\t) = ($^(t)<l>^(t)) , (4.16) 

where the latter expression is defined for n ^ 1. At the one- loop level there are no contribu- 
tions from the q ^ noise to (|4.15| ) and ( [4.16| ). Hence, we have simply subsituted 0o — »• $ 
in these expressions (see ( p.l2| ) and Appendix A). The variances (|4.15| ) and ( [4.16|) are nu- 
merically somewhat easier accessible, yet display a qualitatively similar time dependence as 
^{t). Restricting ourselves again to the critical case r = 0, we have the scaling forms 

m^'\t,L,mi)^L^-^T{t,t,) (4.17) 

and 

m'-l\t,L,mi)^L^-^T^{t,U) . (4.18) 

In order to gain insight into the asymptotic behavior of the scaling function T, we can 
use the fact that the bulk quantity mP'\t,oo,mi) must vary ~ as t — s> oo with m,- 

fixed, where {2 — rfj/z = 1 on the one-loop level. Accordingly, we should find T (t, tj) ~ t 
in the regime t, -C t -C 1. The question then is whether there is a distinct (anomalous) t- 
dependence in the short-time bulk regime t <^ tj -C 1. We claim that this is not the case. The 
way to see this is to make a short-time expansion. As discussed in Refs. f^J^, the leading 
operator contributing to the short-time expansion of 0(a;, t) (in correlation and response 
functions whose other time arguments all stay away from zero) is dt4>{x, t = 0) because 
of the Dirichlet initial condition [§,0. In our case the time argument t of both external 
points approaches zero. However, graphs contributing to quantities such as {dt4>\t=o dt(p\t=o) 
necessarily vanish since the causal structure of the theory leaves only intervals of width 
zero for the integration over internal time arguments. Hence the contribution from the one 
operator in the short-time expansion of m*^^^ should give us the limiting form of T in the 
above short-time bulk regime. This means that the algebraic t-dependence T ~ t should 
remain the same irrespective of whether tj <^ t or tj ^ t. 

Detailed analytical information about the variance m\^{t), defined in ( 4.16|) , is available 



in the limit n ^ oo. From the results of Ref. it is straightforward to show that in the 
bulk {d = 3) and for oo 



(21 1 + €0171^1/3 

m^l^=cit J ' ' , 4.19 
1 + C2 m t 



where ci, C2 are nonuniversal constants. For the scaling function T_L(t, tj), defined in ( [4.18[) , 
this means that 

Tx(t < 1, ti = rrii-^/"' < 1) ^ t('-^)/^ ^^^l^l'i ^^ 20) 

1 + m2 1 ^ ' 

with {2 — rf)/ z = 1 and z/xi = 2 in the large-n limit. Thus, in accord with the general 
discussion in the foregoing paragraph, we obtain ~ t for t <^ ti -C 1 in this special case. 
For t ~ tj, on the other hand, a crossover sets in to a limiting behavior of the form T_l ^ t/3 
for tj <^ t <^ 1. In the finite system this crossover should be observable as long as tj < 1. 
For our numerical solutions, we expect a qualitatively similar behavior, in particular, for 
n = 20. 
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Fig. 5.: Scaling function T for the 
variance defined in (4.17) for n = 1 
and the same values of rrij used in 
Fig. 2, with rrij increasing from top 
to bottom. The smaller figure in- 
serted shows the data in double- 
logarithmic representation. The 
dashed lines indicate the bulk power 
laws ~ t for rrij — > (upper line) and 
rrij — > oo (lower line). 



1 



3 4 5 




Fig. 6.: Scaling function T defined 
in ( 4.17| ) for n = 20 and the same 
values of rrij used in Fig. 3. rrij is 
increasing from top to bottom. 



Fig. 7.: Scaling function T_|_ de- 
fined in ( 4.1 j ) for n = 20. rrij is 
increasing from top to bottom. Ad- 
ditionally to the values for rrij used 
in Fig. 3, one set of data for rrij = 60 
has been added. 
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The numerical data for T and are displayed in Fig. 5 for n = 1 and in Figs. 6 and 
7 for n = 20. The results are in conformity with the above discussion and, in the case of 
n = 20, with the exact large-n result. In the log-log plots inserted the asymptotic bulk cases 
~ t are indicated by the dashed line. The upper and lower dashed lines respresent the limits 
rrij ^ and rrij —>■ oo, respectively. In Figs. 5 and 6, these lines have been fitted to the data. 
In Fig. 7, however, the amplitude of the lower line is one third of that of the upper line, 
which demonstrates that the data for n = 20 are described quite well by the exact solution 



V. CONCLUDING REMARKS 

Although critical relaxation of systems with nonconserved order parameter from a given 
initial state is a long-studied problem 0,^-0 , a full understanding of this phenomenon and 
a theoretical description covering all time regimes with universal behavior has only emerged 
recently for bulk systems In particular, serious efforts to incorporate correctly 



the intial condition and to investigate its long-time traces have only been undertaken in the 
past few years. 

In the present work we have used a combination of field theory and numerical methods 
to study critical relaxation in a finite system with cube geometry and periodic boundary 
conditions. Since we have taken care to build in the initial condition, our analysis — unlike 
previous approaches — has encompassed a proper description of the short-time anomalous 
behavior of bulk systems found by Janssen et al. on the one hand, and on the other, has 
enabled us to investigate quite generally the dependence of the relaxational behavior on the 
initial magnetization rrii in the various (macroscopic) time regimes in a systematic fashion. 
Our main results are the order parameter scaling functions a(t, U) and the scaled amplitudes 
(Too(ti) depicted in Figs. 2 and 3, and in Figs. 4, respectively. Our results for n = 20 (Fig. 3 
and 4) are in very good agreement with the exact n — oo results of our previous work ||10|| . 



Our results for n = 1 display the same qualitative behavior as those of Ref. [|1^ obtained 
by Monte Carlo simulations of the three-dimensional Glauber model. The agreement with 
the latter data cannot be quantitative, of course, simply because the one-loop values of the 
critical exponents pertaining to our approximate calculation differ for e = 1 from the d = 3 
Monte Carlo values. 

Our findings corroborate that the initial condition under favorable conditions leaves 
traces on long time scales. In particular, we have been able to confirm that the amplitude of 
the order parameter m{t) in the regime of linear relaxational behavior exhibits a nontrivial 



dependence on the initial magnetization of the kind found in Ref. [jTO[- Such long- 
time signatures of the initial nonequilibrium state are reminiscent of aging effects (see, e.g., 
Ref. ||2^ for a recent list of references). Aging effects — i.e., the persistence of out-of- 



equilibrium effects after very long times — were originally thought to occur only in systems 



with frustration. As pointed out recently by Cugliandolo et al. |2^, such effects appear 
already at the level of the simplest nonequilibrium systems such as the random walk or the 
Langevin equation for a free scalar field. (Part of the findings of these authors for the free 
scalar field are implicitly contained already in the papers of Janssen et al. |^ and Janssen 
n, which — surprisingly enough — were apparently overlooked.) 
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Aside from being interesting from a conceptional point of view, the improved under- 
standing of long-time traces of the initial condition in critical relaxation from non-equlibirum 
states that has been gained may have potential benefits for numerical investigations of critical 
behavior. As discussed recently by Li et al. |3^, Monte Carlo simulations of the short-time 
anomalous behavior seems to lend itself well to a determination of bulk critical exponents. 

We have considered in this paper only the case of periodic boundary conditions. It 
would be interesting to extend the present approach to systems with free surfaces. In such 
systems, the short-time anomalies become spatially non-homogenous. So far there exist only 
two studies of short-time anomalies in semi-infinite systems [^,^, which both use the e 
expansion and Monte Carlo methods. None of these studied finite-size effects. 

It would also be interesting to extend our analysis to the case of quenches beneath the 
critical temperature. To cope with this case, our techniques would have to be combined 
with the modified handling of the zero mode for T < suggested by Esser et al. @]. 
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APPENDIX A: LOOP CORRECTIONS TO THE STOCHASTIC EQUATION FOR 

THE ZERO-MOMENTUM MODE 



In this appendix we consider the loop corrections to the stochastic equation ( |2.13| ) for 
the zero-momentum part Upon insertion of the decomposition ( |2.6|) , the stochastic 



equation (|2.1| ) can be split into its q = and q parts 

+Ko"([0o,V^],t)=S"(t) 

and 



(Al) 



(A2) 



where Kq and K'^ denote in abbreviated form the deterministic "force" for the q = and 
q modes, respectively. 

In order to avarage over the q ^ noise i), we introduce the generating functional 



Z[E,J] 



exp 



dtJ{t)(t)Q{t) 



(A3) 



This can be represented as a functional integral. A standard procedure (see, e.g., Refs. 
|J2§) yields 



J] oc / V(t)oV<PoV^V^lj j V'dPl^] X 

lo {E-<i)o-Ko)+i^{^-^-K^)+J<j)o 



exp <5 I dt 



(A4) 



where P['d] represents the Gaussian weight for the q noise. A spatial integration over 
the ip part is implied, and the vector indices are omitted. Contributions resulting from the 
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functional determinant det{S(/S(j)) have been omitted since they vanish if we use a prepoint 
discretization. 

The integration over ^ is elementary, giving 



with 



Z[E,J] 



V(j)oV(poViljV'iljexp{-Jes+ 1^ dt 



Jes = rrii 



0o(O) + dt {00 (00 + Ko) + ^ + K^) - Ao^'} . 



(A5) 



(A6) 



From ( |A5|) we see that S plays the role of a source for 0o in this representation. Thus 
Z[E, J] may be regarded as a generating functional for expectation values of products of the 
operators 0o and 0o- Let us introduce its Legendre transform 



r[l>,$] = -lnZ[S,J] 



where it is understood that the solutions to 



dt E<l> + M 



~, , 6\nZ\E,J] , 5\nZ[E,J] 



6E{t) 



6J{t) 



(A7) 



(A8) 



are substituted for S and J on the right-hand side of ( |A7|) . In terms of r[$, $], the generating 
functional of the one-particle irreducible functions pertaining to the cumulants generated by 
lnZ[S, J], the desired stochastic equation ( |2.13| ) simply becomes 



i{t) 



(A9) 



$=0 



Loop corrections to the "classical" stochastic equation for ^(t) (corresponding to the 
substitution F — > J^s in ( |A9| ) and given by ( p.l3| ) with the terms involving C""^ and C"^'^ 
omitted) can now be systematically computed in a standard fashion. To this end we insert 
(p.l2| ) and the analogous decomposition for (po into J'es- From the terms quadratic in 

ip, 6(j), 6(j) of the action ( |X6| ) we obtain the free response and correlation propagators. 
Their spatial Fourier transforms read, in matrix notation. 



free 



G(q;t,t') ^ (0,(t)0_,(t')>7^ = ^(t-t') exp{-Ao^, dt" \^^q^ + + f M{t")^'j 



(AlO) 



and 



mm{t,t') 



C{q- 1, t') ^ {i:,m-,it'))'r = 2Ao / dt" G{q, t, t") ■ G{q, t', t") , (All) 



where Al was defined in ( p.32|) . While ( |A10|) also holds for q = 0, the free correlation 
propagator {5(f){t)5(j){t')y^'^'^ vanishes identically. The latter is due to the absence of a term 
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quadratic in 0o in iTefr, which in turn is a direct consequence of the fact that we averaged 
only over the q ^ part d or the noise, but not over its q = part S. 

Up to two-loop order the functions C'^^{t) and C'^^'^{t) in ( |2.13| ) are given by 



0.6 



and 



+ 0(3-loops) 



(A12) 




0(3-loops) 



(A13) 



where the full lines with and without an arrow denote the free response and correlation 
propagators (|A10| ) and ( |A11| ), respectively, while the broken line signifies 
The explicit expressions for the graphs are 



O 



(AM) 




XogoL-'' j dt' t, t') G^^\q; t, t') S'^'f''^^ C^\t') 





(A15) 



and 



t = 3Ao(7o^"'' 5"'^'^''' / dt' it') C'""' {q- 1, t') C^^' {q'; t, t') G^^' {q"; t, t') 



q + q' + q"=0 



(A16) 



From expression ( |A14| ) the one-loop result ( |3.29| ) for C used in the main text follows in 
a straightforward manner. 
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APPENDIX B: LAURENT EXPANSION OF REQUIRED INTEGRALS 

Equations ( |3.34| ) and ( |3.35D involve the integrals 



ds 



AVks) 



-pKS 



(Bl) 



with K = Sh'^XqL^^, p = L^tq/Atc'^, and p = 0,1. Here we will compute the Laurent 
expansion of these to 0(e°). 

From A'^{ks) we add and subtract its limiting form (Ks/vr)"^^'^ for s — 0. This gives 

= [n'^/^JpiKt, p) + Fp{Kt, p) + 0(e)] , (B2) 

where Fp is the function defined in ( |3.4(j| ), while 

nv,p)= r dxxP'^^'/^ e-""" = T{p - I + e/2) yP-'^'/^ -f*{p ~ 1 + e/2,vp) . (B3) 
Jo 

Here 7*(a, x) = x~"7(a, x)/r(a) is the analytic function related to the incomplete gamma 
function 7(0, x) The poles arise solely from Jp. Therefore we have expanded the second 
(regular) term inside the brackets. 

The Laurent expansion of the function Jp can be obtained in a straightforward fashion, 
utilizing the recursion relation 



7* -1 



x'y 



X + 



exp(— xj 



familiar properties of r(a), and the expansion 

e 



1 y-^.^ 



1 + 2 



Ei(-x) -Inx +C(e' 



(B4) 



(B5) 



where the exponential-integral function Ei(a;) of Ref. [^] used here is related to the function 
El of Ref. g via Ei(-a;) = -Ei(a;). 
One finds 



Jo{v,p) = p 



-2 
e 



1 + In p — Ei(— t;p) — 



exp(— t;p) 
vp 



+ 0{e) 



and 



Mv,p) = - + Ei{-vp) -\np-CE + 0{e) . 



(B6) 



(B7) 



Substitution of these results into (M^ leads to 



Io{t) = X,'KdLUL'ro 



-2 , 2Aot exp(-2Aorot) 



and 



hit) 



+ 0{e) 



(B8) 



TT 



d/2 



- + B{2\oTot) + In L"^ To- Ce 



(B9) 
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